Does Low Frequency X-ray QPO Behavior in GRS 1915+105 
Influence Subsequent X-ray and Infrared Evolution? 



Valerie J. Mikles 

Department of Astronomy, University of Florida, Gainesville, FL 32611 

mikles@astro . uf 1 . edu 
Stephen S. Eikenberry 
Department of Astronomy, University of Florida, Gainesville, FL 32611 

eikenberry@astro . uf 1 . edu 
and 

David M. Rothstein 
Department of Astronomy, Cornell University, Ithaca, NY 14853 
droth@astro . Cornell . edu 

ABSTRACT 

Using observations with the Rossi X-ray Timing Explorer, we examine the 
behavior of 2 — 10 Hz quasi-periodic oscillations (QPOs) during spectrally-hard 
dips in the x-ray light curve of GRS 1915+105 that are accompanied by infrared 
flares. Of the twelve light-curves examined, nine are /3-class and three are a-class 
following the scheme of Belloni et al. (2000). In most cases, the QPO frequency 
is most strongly correlated to the power law flux, which partially contradicts 
some earlier claims that the strongest correlation is between QPO frequency 
and blackbody flux. Seven /3-class curves are highly correlated to blackbody 
features. In several cases, the QPO evolution appears to decouple from the 
spectral evolution. We find that /3-class light-curves with strong correlations can 
be distinguished from those without by their "trigger spike" morphology. We 
also show that the origin and strength of the subsequent infrared flare may be 
causally linked to the variations in QPO frequency evolution and not solely tied 
to the onset of soft x-ray flaring behavior. We divide the twelve a- and /5-class 
light-curves into three groups based on the evolution of the QPO, the morphology 
of the trigger spike, and the infrared flare strength. An apparent crossover case 
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leads us to conclude that these groups are not unique modes but represent part 
of a continuum of accretion behaviors. We believe the QPO behavior at the 
initiation of the hard dip can ultimately be used to determine the terminating 
x-ray behavior, and the following infrared flaring behavior. 

Subject headings: accretion, accretion disks- black hole physics- stars: individual 
(GRS 1915+105)- stars: oscillations 

1. Introduction 

Discovered in 1992 by Castro-Tirado et al. (1992), GRS 1915+105 is an x-ray transient 
that continues to intrigue us with its unique array of variability on many wavelengths and 
timescales. Dubbed a microquasar because of its apparent superluminal jets (Mirabel & 
Rodriguez 1994), GRS 1915+105 is a black hole candidate and x-ray binary. Although it 
suffers ~ 20 — 30 magnitudes of extinction at visible wavelengths, GRS 1915+105 has shown 
great activity in the radio, infrared, and x-ray regimes. Of particular interest in this study 
are x-ray light-curves showing spectrally-hard dips on ~ 30 minute timescales. 

It is well-established that the spectrally-hard dips that frequently appear in the x-ray 
light-curves of GRS 1915+105 are associated with infrared and radio flares (Eikenberry et al. 
1998; Mirabel et al. 1998; Fender & Pooley 1998; Klein- Wolt et al. 2002). These spectrally- 
hard dips are also associated with 2 — 10 Hz variable quasi-periodic oscillations (QPOs). 
Belloni et al. (2000) defines a total of twelve x-ray light-curve classes for GRS 1915+105, 
distinguishable by general appearance, count rate, and x-ray color (related to hardness). We 
study the a- and /3-class light-curves in the 2 — 25 keV range. Both have extended spectrally- 
hard dips but the length of the a-class hard dip is nearly twice the length of the /5-class. 
The a-class light-curves have a dip of ~ 1200 seconds followed by strong x-ray oscillations 
marking the end of the hard state. The /5-class light-curves are more complex. The dip 
lasts 500 — 700 seconds and its ending is marked by a spectrally soft spike followed by a 
nearly monotonic rise. After reaching a peak flux, the x-rays begin rapid, large-amplitude 
oscillations. In both the a- and /5-classes, simultaneous x-ray and infrared observations show 
infrared flares rising as the hard x-ray dip ends (Eikenberry et al. 1998; Mirabel et al. 1998; 
Rothstein et al. 2005). In the /3-class, the infrared flare begins near in time to the soft spike, 
called the "trigger spike". Eikenberry et al. (1998) observed a one-to-one correspondence of 
x-ray dips to infrared flares at 30 minute intervals. 

While the x-rays are rapidly oscillating, the infrared flare peaks and decays without the 
rapid, large amplitude variations seen in the x-ray. Flares associated with our observations 
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are defined by Eikenberry et al. (2000) as class C and class B, ranging from 30 — 200 mJy 
when de-reddened by 3.3 magnitudes. Individual Class C infrared flares are smaller (5 — 10 
mJy) and have been observed individually associated with isolated x-ray flares in a dip/flare 
cycle (Eikenberry et al. 2000). When several class C infrared flares occur in rapid succession, 
they can appear as a larger flare. Rothstein et al. (2005) shows this to be the case for the 
a-class light-curves in our sample, where the associated infrared flares range from 10 — 30 
mJy. Rothstein et al. (2005) showed that if each soft x-ray flare were associated with a 5 — 10 
mJy (class C) sub-flare then the duration and strength of the overall infrared flare would 
be explained. By that analysis, the predicted collection of infrared sub-flares associated 
with a /3-class light-curve would contribute only a small fraction of the observed 60 — 200 
mJy infrared flux. This small contribution was observed by Eikenberry et al. (2000) as an 
"infrared excess" . Although the period of the rise and decay of the primary class B infrared 
flare is not coupled to the period of the x-ray oscillations, the duration of the infrared excess 
is (Eikenberry et al. 2000; Rothstein et al. 2005). The difference between the dip/flare cycles 
associated with a- vs. /3-class curves becomes the presence of a large primary infrared flare 
which is uniquely associated with /5-class x-ray light-curves. The presence of the primary 
flare is widely associated with the presence of the x-ray "trigger spike" . It is believed that 
both are linked to the underlying cause of larger plasma ejections. 

Several authors have shown that infrared flares tend to be followed by radio flares 
of similar morphology (Fender & Pooley 1998; Mirabel et al. 1998). The sequence of a 
spectrally-hard dip, an infrared flare, then a radio flare, is generally associated with a plasma 
ejection from the source. The ejection is observed as the x-rays transition from a spectrally- 
hard to a spectrally-soft state. During the hard dip that precedes an ejection event, a 
2 — 10 Hz QPO is always observed (Belloni et al. 2000). While the QPO is present, the 
disk x-ray emission is greatly reduced in the 2 — 25 keV range and the x-ray luminosity is 
dominated by power law flux. Several groups have found the QPO peak frequency positively 
correlated to both power law and thermal disk components, suggesting that the QPO arises 
in the same location as the related emission or is causally related to it (Muno et al. 1999; 
Feroci et al. 1999). Markwardt et al. (1999) observed that at frequencies above 4 Hz, the 
QPO is most strongly correlated to the thermal disk component, specifically the blackbody 
disk flux. At lower frequencies, the QPO shows a broad association with the power law 
flux. Muno et al. (1999) also show the QPO is strongly associated with disk temperature 
during the hard dip. After the trigger spike, the power density spectrum becomes smooth 
(Markwardt et al. 1999). 

In this paper, we examine the behavior of the QPO during the hard dip and the relation- 
ship of its behavior to subsequent infrared flares. In Section 2, we discuss our observations 
and analysis technique. We calculate and discuss the correlation of the QPO peak frequency 
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to spectral features in Section 2.1 and comment on the time evolution of the QPO during 
the dip in Section 2.2. In Section 2.3, we examine the morphology of the trigger spike in the 
/3-class light-curve and the indication of a continuum of behaviors from this class to a-class, 
which does not show a trigger spike. In Section 2.4, we compare infrared flaring behavior 
between the events. We discuss the significance of these results in Section 3 and in Section 
4 we summarize our conclusions. 



2. Observations and Analysis 

We use several Rossi X-ray Timing Explorer (RXTE) observations of GRS 1915+105 
taken on 14 - 15 August 1997, 9 September 1997, and 27 - 28 July 2002 (see Table 1). 
We extract Proportional Counter Array (PCA) Standard-1 light-curves using FTOOLS v5.3 
and identify regions showing a hard x-ray dip. We examine 12 regions in this work; nine of 
them are /3-class and three are a-class (Belloni et al. 2000). They are numbered sequentially 
f3-l through (3-9 and a-10 through a-12. For the twelve regions, we extract binned mode 
8-millisecond light-curves in the 2 — 13 keV range and 4-second resolution binned and event 
x-ray spectra in the 2 — 25 keV range. Using XSPEC vl 1.3.0, we fit the spectra with a 
combination of absorbed multi-temperature disk blackbody and power law models. Only fits 
with xl < 2 are used in analyses. 

Using a Fast Fourier Transform, we calculate the power density spectrum (PDS) from 
the binned mode 8-millisecond light-curve. We fine-bin the PDS using Fourier interpolation 
and track the peak frequency at 4-second resolution (see e.g. Ransom et al. 2002). Fine 
binning allows us to calculate a higher-resolution Fourier response by interpolating responses 
at non-integer frequencies. In Figure 1, we show a gray-scaled fine-binned PDS and overplot 
a one-second resolution light-curve. 

We determine the QPO frequency by fitting a Moffat function to the PDS in the 2 — 10 Hz 
frequency range. The Moffat function is a Lorentzian modified with a variable power law 
index (Moffat, A. F. J. 1969). The QPO is considered detected if it has a quality factor 
Q = u/FWHM > 2, where v is the centroid frequency of the Lorentzian and FWHM is 
the full-width at half max. In a-class curves, we require the QPO frequency to be between 
2.6 and 10 Hz to avoid contamination by low frequency noise. In /3-class curves, we allow 
detections in the full 2 — 10 Hz range. However, also visible in the PDS is a low-frequency 
noise component, the maximum frequency of which tends to rise above 2 Hz at the beginning 
and end of the dip. To avoid spurious detections from this component, we require that each 
detected QPO does not vary too sharply from the QPOs at surrounding times. For points 
at times t > 350 seconds, we take an average of the prior 10 detected frequencies, v avg . The 
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next detection is required to be greater than 75% of v avg . For t < 150 seconds, QPOs are 
treated similarly, though u avg is determined by using points following the reference point, as 
opposed to those preceeding it. Although omitting points as non-detections can leave part 
of the QPO evolution under-sampled, we believe we have sufficient representation from the 
regions to ensure a sound qualitative result and a reasonable quantitative one. In addition 
this method yields a consistent and repeatable QPO detection. It is possible that adjusting 
the range in this manner causes us to omit real, unresolved oscillations in the QPO, but visual 
inspection shows that these points are not part of the primary U-shaped QPO feature we 
wish to focus on. This shape is apparent in the PDS shown in Figure 1. In Figure 2, we plot 
the detected QPOs as filled circles against a 1-second resolution light-curve. For illustrative 
purposes, we show power peaks in the PDS that do not meet the detection criteria as open 
circles. 



2.1. QPO Frequency - Correlation to Spectral Features 

For each hard dip, we calculate the Linear Pearson Correlation Coefficient, r, between 
the QPO frequency and various spectral features. We define events with |r| > 0.70 as highly 
correlated. This value is chosen as it means that at least half the variance of the spectral 
feature can be accounted for by the variance of the QPO frequency. Values of |r| between 0.4 
and 0.7, we discuss as believable, but generally disregard when more statistically significant 
correlations are present. 

We list the correlation coefficients for the twelve dips in Table 2. In all of the cases, 
we observe a highly significant correlation between the QPO frequency and the total flux 
which is generally stronger than those to individual blackbody or power law features. When 
considering model-specific spectral features, eleven out of twelve show strong correlations 
to power law flux, while only six out of those eleven show strong correlations to blackbody 
features. In nine out of those eleven cases, the correlation to the power law flux is stronger 
than the correlation to any blackbody feature. This gives the apparent result that the 
QPO frequency is more fundamentally tied to the power law component than the blackbody 
component. 

Figure 3 shows a scatter plot of the QPO frequency versus the power law flux for the 
twelve cases. The points marked with triangles are detections within the first 100 seconds 
of the dip and the squares are from the last 100 seconds of the dip. The line is based on a 
fourth order fit to the time evolution of the QPO frequency and the power law flux, tracing 
the approximate path of the evolution. The open circles are power peaks between 2 and 
10 Hz where a QPO was not detected (see above). These points are generally not a part of 
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the observed trends. 

In (3-1 through (3-5 (August 1997) the evolution of the QPO-power law flux relation is 
generally tighter and the correlations stronger (r > 0.80). The a- 10 through a-12 curves also 
show a tight correlation, though the QPO frequency varies over a smaller range of frequencies. 
Like the a-class curves, we only have detections for (3-8 over a small range of frequencies, 
during which we see a tight correlation to the power law flux. In the other three cases, (3-Q, 
(3-7, and (3-9, we notice a hysteresis effect, observed as a separation between the entrance 
(triangles) and exit (squares) points of the dip. In the (3-7 case, the hysteresis is extreme and 
lowers the correlation coefficient drastically despite the initially linear correlation at the dip 
entrance. The squares which trace the last 100 seconds of the dip show the QPO frequency 
rising sharply while the power law flux is relatively constant, suggesting a decoupling of the 
features. 

In (3-1 through (3-7, the correlation between the QPO frequency and the power law flux 
is accompanied by a similarly strong correlation to the blackbody flux and blackbody disk 
temperature. In Figure 4, we see a strong correlation to the blackbody flux is most prominent 
above 4 Hz, a point also observed by Markwardt et al. (1999). However, the correlation does 
not last over the broad part of the hard dip where the QPO can change significantly while 
at relatively constant blackbody flux. In (3-8 and (3-9, where the blackbody flux appears 
uncorrelated, there are fewer QPO detections above 4 Hz due to the rise in low frequency 
noise at the dip exit. 

We show the QPO frequency - blackbody temperature scatter plots in Figure 5. Clearly 
the QPO in the a-class curves is not believably correlated to the blackbody temperature. 
In (3-1 through (3-5 and (3-8, the correlations to the blackbody temperature follow a steady 
trend, but suffer a slightly wider dispersion in their evolution than the power law flux. In 
(3-6, (3-7, and (3-9, we see a similar hysteresis to that observed in the power law flux. To 
untangle the possible interplay of the blackbody temperature and power law flux, we apply 
a partial correlation analysis. The partial correlation coefficients are listed in Table 3. In 
the table, we calculate coefficients for four scenarios: 

1. The QPO frequency — power law flux correlation removing the effect of blackbody 
flux. 

2. The QPO frequency — power law flux correlation removing the effect of blackbody 
temperature. 



3. The QPO frequency — blackbody flux correlation removing the effect of power law 
flux. 
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4. The QPO frequency — blackbody temperature correlation removing the effect of power 
law flux. 

In the first scenario, we see that with the blackbody flux removed, the correlation to 
power law flux is still strong. This is expected because the blackbody flux poorly explains 
variances at frequencies below 4 Hz (see Fig. 4 and Table 2). In the second case, we 
find the removal of blackbody temperature has a more significant effect. In most cases, 
the QPO frequency — power law flux correlation is weak, though believable (r > 0.4). 
This means that after removing variations in the QPO frequency and power law flux that 
can be explained by variations in blackbody temperature, variations in the QPO frequency 
remain that can be at least partially explained by variations in the power law flux. High 
significance is seen in a-class cases, which is expected due to their weaker dependence on 
blackbody temperature. In the third case, we test the correlation to the blackbody flux 
after removing the power law flux. These correlations are believable and occasionally strong, 
suggesting that a combination of the blackbody and power law sources is required to explain 
the QPO frequency. In the final case, the removal of power law flux from the QPO frequency 
- blackbody temperature correlation, we see that partial correlation coefficient drops below 
significance in half of the cases. This means that once variations in the QPO frequency that 
can be explained by variations in the power law flux are removed, no variation remains that 
can be explained by blackbody temperature. 

In summary, the QPO frequency is often correlated to the power law flux. For /3-class 
light curves, when this correlation is strongest, we tend to find a correlation to the blackbody 
flux and blackbody temperature as well. In cases where there is a slightly weaker correlation 
to the power law flux, the correlation to the blackbody features is less predictable and a 
hysteresis effect is visible in the power law flux and blackbody temperature relations. In 
contrast, a-class QPOs tend to have a strong correlation to the power law flux and weak or 
non-existent correlations to the blackbody features. Based on this evidence, we believe the 
correlation between the QPO frequency and the power law flux is more fundamental. 

2.2. QPO Time Behavior 

Noting the hysteresis in the QPO frequency versus power law flux distribution in several 
of the /3-class light curves (see Fig. 3), it is likely that the spectral evolution is intimately 
tied to the QPO time evolution. The QPO evolution at the beginning of all the hard dips 
is similar, each beginning with the sudden appearance of a QPO at 6 — 10 Hz. This QPO 
will smoothly drop to between 2 — 3 Hz within seconds of the initiation of the hard dip. 
After the x-ray trigger spike, the primary U-shaped QPO feature will disappear, replaced 
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by occasional power peaks. In Table 4 we list the minimum QPO frequency during the dip 
and the time spent near that frequency. Based on the observable variations, we then divide 
the light-curves into three broad groups: 



• Group 1: Figure 6 shows an example of the /3-class light-curves in Group 1. A series 
of x-ray oscillations calm into a low, spectrally-hard dip within about 100 seconds. 
During this time, a QPO arises at ~ 6 — 8 Hz. Following the intensity drop of the 
light-curve and more specifically the power law flux, the QPO falls steadily to ~ 2 Hz. 
It remains at this frequency for over 150 seconds (see Table 4), after which the power 
law flux and QPO frequency begin a slow rise. The U-shaped QPO vanishes after the 
x-ray trigger spike. The total length of the dip is on the order of 600 seconds, and ~ 
30% of that time is spent at the minimum frequency. 

• Group 2: This group is also composed of /5-class light-curves with similar spectral 
behavior to Group 1 (see Fig. 7). However, the QPO behavior in this group is some- 
what different. While in Group 1, the QPO falls off to ~ 2 Hz and lingers, in Group 2 
the QPO immediately starts to rise again. In Table 4, we show that while the Group 1 
events remain near the minimum frequency for > 150 seconds, the Group 2 events re- 
main for < 100 seconds (about 15% of the dip length). The following rise in frequency 
is the start of the hysteresis in the QPO frequency — power law flux scatter plot (Fig. 
3) and likely indicates that the QPO and power law flux have decoupled. In addition, 
these cases see the rise of a low frequency noise component above 2 Hz as the QPO 
weakens in amplitude and rises rapidly in frequency. While these low frequency points 
are excluded from correlation analysis as being associated with low frequency noise, it 
is possible that they represent an increase in rapid, unresolved QPO oscillations. The 
difference in behaviors in the first two groups is most remarkable because their x-ray 
light-curves and spectral behaviors are so similar. 

• Group 3: The final group contains the a-class light-curves represented in Figure 8. 
In these events, the hard dip is surrounded by x-ray oscillations but no independent 
terminal spike is observed. The total length of the dip is ~ 1200 seconds and the QPO 
disappears when the dip ends. This case is similar in shape to Group 1 QPO evolution, 
though twice as long. It differs in that on entering the dip the QPO frequency levels 
off at ~ 3 Hz. Overall, the frequency varies over a smaller range than seen in the other 
two groups. 
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2.3. Differing Trigger Spike Morphology in /3-class Light-Curves 

The separation of the /3-class curves into two groups of QPO evolution inspires the 
search for other differences in the light-curve behavior. We focus here on the trigger spike 
which appears at the end of the hard dip and coincides with the start of the infrared flare and 
note that the intensity and morphology of the spikes vary between Group 1 and Group 2 
events. Figure 9 shows a one-second resolution light-curve for each of the nine observed 
spikes in a time range of 55 seconds before and after the maximum. 

The light-curves are fit with a double Gaussian plus a polynomial of the form: 

/ = A e- Z * + B e~ z B + C + C x t 
where Za = (t — t a )/ct a and z B = (t — r B )/a B . 

The second Gaussian term fits a small rise preceeding the x-ray spike - - a feature 
more distinct in the (3-3 through (3-7 light-curves. This second Gaussian is treated with the 
polynomial as part of the background emission. A normalized peak amplitude for the spike 
is calculated as A = Aq/A^ where A^ g is the polynomial and second Gaussian evaluated at 
the primary Gaussian peak time, t = ta- It should be noted that the apparent symmetry 
of the spike is tied to the time resolution. The one-second time resolution allows a fair 
sampling of data points for the range while smoothing the rapid x-ray variability seen in 
an eight-millisecond curve. Broader time intervals will leave the spike under-sampled and 
reduce the accuracy of the fit. 

The relevant numerical fits are listed in Table 5. It is interesting to note that the 
integrated count rate (/ iTit ) over the full-width-half-max of the peak is similar for all data 
sets, but in most other aspects of the fit, the August (Group 1) and September (Group 2) 
data have different properties. 

• The Group 1 data have a higher normalized amplitude: A Grpl > A Grp2 . 

• The Group 1 data are more symmetric while the Group 2 data shows a sharp cut-off 
after the peak flux (see Fig. 9). 

• The Group 1 data spikes are narrower: cr Grpl < <7 G A rv<1 . 

• The underlying slope of the light-curve is positive in the Group 1 data and flat or 
negative in the Group 2 data: C Gtpl > while C Gtp2 < 0. 

While all Group 2 events show a negative underlying light-curve slope, it is interesting 
that two of them ((3-6 and (3-7) have a nearly flat slope while the other two ((3-8 and (3- 
9) have a decidedly negative slope. From Figure 3 we note that the two with nearly flat 
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slopes have a slightly more visible hysteresis because they have more QPO detections above 
4 Hz. The low frequency noise component does not rise as strongly above 2 Hz and they 
show a stronger correlation to blackbody temperature. In addition, the terminating spikes 
of (3-8 and (3-9 look much more disturbed than those of other light-curves (Fig. 9). The 
differences in trigger spike morphology suggest that the variation in QPO behavior is not 
an artifact of the detection method. The (3-6 event is particularly interesting because it 
can be considered a crossover between Group 1 and Group 2 events. The QPO frequency 
- power law flux relation in (3-6, shows a significant correlation, r = 0.78, despite the 
apparent hysteresis. From Figure 3, we see that this event seems to be a bridge between the 
sharp linear correlations of Group 1 and the divergent shapes of Group 2. We classify it as 
Group 2 because the QPO frequency clearly deviates from the initial regression and because 
of its trigger spike morphology. The trigger spike of (3-6 is wider and more asymmetric 
than Group 1 events. In addition, the flat underlying light-curve of (3-6 suggests that it 
is more appropriately associated with Group 2 than Group 1 events. This being said, we 
acknowledge that the groups are not absolute, but likely part of a continuum of behaviors. 



2.4. Associated Infrared Flaring 

Simultaneous infrared coverage is available for four out of five of our Group 1 data sets. 
Eikenberry et al. (1998) show that these dip/spike pairs are usually followed by large infrared 
flares. The events observed ranged from ~ 60 to 200 mJy Although the rise and fall of 
the flare does not correspond to the period of x-ray oscillation, a weak infrared excess which 
lasts throughout the period of the x-ray oscillations is observed. Eikenberry et al. (2000) 
and Rothstein et al. (2005) explain this excess as the superposition of many faint infrared 
flares each on the order of ~ 10 mJy. The dominant infrared flare is associated with the 
x-ray trigger spike. 

Mirabel et al. (1998) observed a single infrared flare event associated with our (3-7 curve 
which reached an amplitude of ~ 30 mJy. Simultaneous infrared coverage is not available 
for the other hard dips in Group 2. In all three a-class light-curves of Group 3, the hard dip 
is followed by a ~ 30 mJy infrared flare (Rothstein et al. 2005). In Mirabel et al. (1998), 
an a-class event is observed to be followed by a flare that peaked at ~ 10 mJy. Rothstein 
et al. (2005) showed that the ~ 30 mJy flares can be explained as a summation of Class C 
sub-flares, each associated with a soft x-ray flare. 
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3. Discussion 

3.1. QPO Correlation with Spectral Features 

Previous studies have shown that the QPO is most strongly tied to the thermal disk 
component (Markwardt et al. 1999; Feroci et al. 1999; Muno et al. 1999). Using the Septem- 
ber 1997 data (our (3-6 through (3-9), Markwardt et al. (1999) point out that the correlation 
to disk flux is strongest when the QPO frequency is above 4 Hz. While this is true, the 
QPO frequency is in this range less than 25% of the time, and mostly falls into this range 
when entering or exiting the hard dip. During the course of the hard dip, the QPO will 
change significantly while the blackbody flux remains relatively constant. Markwardt et al. 
(1999) also say that at lower frequencies, there is an apparent broad correlation with power 
law flux. We confirm that this correlation may exist, and show that it is most apparent 
at lower frequencies. Because of a hysteresis effect, the deviation in the QPO frequency - 
power law flux relationship is more apparent at high frequencies. We suggest that an initially 
tight correlation is broken as the QPO begins to rise in the latter half of the dip. Two out 
of four of the observed September events are strongly correlated to power law flux and a 
different two out of four are correlated to blackbody temperature. 

In contrast, the August 1997 (our (3-1 through (3-5) events show a strong correlation 
with both power law and blackbody features - specifically the power law flux and blackbody 
temperature. A partial correlation analysis shows that if either the power law flux or black- 
body temperature is removed, the correlation to the other is weakened, so it is not likely the 
effects of these components can be untangled. We do, however, argue that the correlation to 
the power law component may be more fundamental, especially since the power law flux is 
more strongly tied to the x-ray emission at this point and the disk component is vanishingly 
small. In addition, the a-class curves show a consistently strong correlation to the power 
law flux and less consistent correlation to blackbody features. 

As mentioned before, the bulk of the QPO change occurs during the entry into and exit 
from the dip. While in the dip, the spectral features remain fairly stable. Thus much of the 
correlation strength is tied to QPO stability during the dip and the relatively short motion 
at the start and end of the dip. This being the case, the rise in low-frequency noise which 
greatly affects QPO detections in (3-8 and (3-9 is also likely to be affecting the correlation 
coefficients. It is reasonable to assume that the number of non-detections late in the dip 
is reducing the number of QPO detections above 4 Hz and thus artificially lowering the 
correlation to blackbody features in these two cases. Short time-scale variations (where the 
QPO suddenly dips and recovers) also tend to create outliers to the correlation because the 
change in frequency is not accompanied by spectral variation. It is uncertain if these jumps 
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are real as the shape of the variation is indeterminant on four-second time scales. These 
outliers are omitted from our correlation analysis. 

3.2. Infrared Flaring Behavior 

When comparing these events, we cannot ignore the fact that all these hard dips are 
followed by an infrared flare, the strength of which is tied to the presence or absence of a 
terminating spike. Furthermore, the shape of the spike is tied to the behavior of the QPO 
during the preceding dip. Therefore, the properties of the QPO appear to be fundamental 
for determining the subsequent jet activity in GRS 1915+105, and the question becomes, 
"If all hard dips start relatively the same (with the appearance of a QPO), why don't all 
hard dips end the same?" What causes the QPO to quickly rise in Group 2 or to have a 
higher minimum frequency in Group 3? Is the energy that would be fed into the QPO and 
subsequent infrared flare for a Group 1 event being leaked out before the dip termination in 
Group 2 events, or is the energy missing from the system altogether? The Group 1 events 
have strong, symmetric trigger spikes and strong ~ 100 mJy infrared flares. The relatively 
weak trigger spike and ~ 30 mJy flare observed in September 1997 may be evidence that a 
weaker trigger spike would be associated with a weaker primary flare. 

On a final note, what we have referred to as the "trigger spike" of the /5-class is so named 
because it coincides with the start of the infrared flare. However, observations of several (3- 
class events by Eikenberry et al. (1998) were unable to conclusively determine whether the 
infrared flare began simultaneously with the trigger spike. Also, in looking at the Mirabel 
et al. (1998) event corresponding to our f3-7, one might believe that the infrared flare starts 
100 — 200 seconds PRIOR to the spike (see Mirabel et al. 1998, their Figure 3). Noting that 
the QPO also significantly weakens compared to the low frequency noise component 100 — 200 
seconds prior to the spike (see Fig. 1) suggests that the origin of the infrared flare may be 
causally linked to the mechanism powering the QPO, and not tied initially to the soft x-ray 
flaring or trigger spike. This suggests that the QPO is tied to a multi-wavelength energy 
release or the formation of mult i- wavelength features. The decoupling of the QPO from x-ray 
spectral features in this case supports this. In all other cases, the disappearance of the QPO 
coincided with the trigger spike or the first x-ray oscillation (and thus the infrared flare), so 
this picture would be consistent. Eikenberry et al. (2000) observed a series of class C infrared 
flares preceding the x-ray soft flares in #-class dip-flare cycles, so this sequence would not be 
entirely unprecedented. 
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3.3. A Cause And Effect Summary 

While the nature of the QPO is still uncertain, it does reflect and may possibly be used 
to predict observable outcomes. In all of these cases, a QPO appears when the x-ray energy 
changes from soft to hard. While the change in hardness occurs very quickly, the flux drops 
slowly for ~ 100 seconds. The QPO, initially around 6 Hz, falls with similar smoothness 
(see Fig. 2). This is where the tracks diverge. We believe that the QPO behavior at the 
divergent point can ultimately be used to predict how the dip will end. Consider the three 
groups we identify, summarized in terms of cause and effect: 



• Group 1: A /3-class light-curve enters a hard dip phase. The QPO falls off to 2 Hz 
and maintains that frequency for > 150 seconds. 

END RESULT: 

- The QPO frequency is tightly correlated to both blackbody and power law spectral 
features and the correlation lasts the length of the dip. 

— The dip lasts 500 — 700 seconds and terminates with a strong, narrow, symmetric 
spike and increasing underlying flux. 

- A strong class B infrared flare of strength ~ 100 mJy follows with infrared excess 
explained by pursuant x-ray oscillations. 

• Group 2: A /3-class light-curve enters a hard dip phase. The QPO falls off to 2 Hz, 
but begins increasing after < 100 seconds. 

END RESULT: 

— The QPO frequency decouples from the spectral features and the strength of the 
correlation is related to the degree of hysteresis. 

— The QPO significantly weakens before the terminating spike. 

- The dip lasts 500 — 700 seconds and terminates with a weak, wide, asymmetric 
spike and flat or decreasing underlying flux. 

- A slightly weaker class B infrared flare of strength ~ 30 mJy follows, possibly 
starting with the weakening of the QPO. 

• Group 3: An a-class light-curve enters a hard dip phase. The QPO falls off to 3 Hz 
and maintains that frequency for several hundred seconds. 

END RESULT: 

— The QPO is correlated to power law features but not to blackbody features. 
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- The dip lasts ~ 1200 seconds and then immediately enters a period of oscillation. 

- A series of class C infrared flares with summed amplitude ~ 10 — 30 mJy follows 
and is associated with the duration of the x-ray oscillations. 

4. Conclusions 

In conclusion, we have identified three groups of hard dips with a range of spectral 
behavior, QPO frequency evolution, trigger spike morphology, and infrared flare strength. 
Each hard dip is associated with a variable low frequency QPO and is followed by an infrared 
flare indicating an ejection event. Similarities in the x-ray light-curve, x-ray hardness, and 
infrared flaring suggest that a similar mechanism is responsible for these behaviors. 

While it is easy to expect different behavior from a- and /3-class light-curves, it is 
surprising to find differences within the /3-class itself. Most interesting is the spectrum 
of QPO time evolution behaviors seen in our small set of observations and the fact that 
only a slight variation is necessary to affect the end result. The correlation of the QPO 
to x-ray spectral features hinges on the time evolution of the QPO. The time evolution is 
also intimately tied to the trigger spike morphology and subsequent infrared flaring. Most 
interesting is the possibility that the infrared flare begins with the disappearance of the 
QPO and is not solely tied to the trigger spike or x-ray flares. It is clear that this study only 
scratches the surface of the spectrum of behaviors exhibited by GRS 1915+105, however, by 
studying these events together, we may better understand the underlying mechanism. 

The authors owe a debt of gratitude to the astute referee for aiding in the refinement 
of this work. V. J. M. and S. S. E. are supported in part by an NSF CAREER Award 
(AST-0328522). V. J. M. is also supported by a University of Florida Alumni Fellowship. 
D. M. R. is supported by a National Science Foundation Graduate Research Fellowship. This 
research has made use of data obtained from the High Energy Astrophysics Science Archive 
Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. 
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Fig. 1. — A gray-scaled, fine-binned power density spectrum for the light-curves. We plot 
the 1-second resolution light-curve over the PDS. Note that, in the (3-8 and (3-9 
strong low-frequency noise component appears ~100 seconds before the spike. 
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Fig. 2. — An overlay of the 1-second resolution x-ray light-curve (line) and fine-binned 4- 
second QPO frequency (circles). Time is in seconds and QPO frequency is in Hz. The open 
circles are low frequency power peaks observed when the QPO is not detected (see Section 
2). They are shown for illustrative purposes. Note that for (3-Q through (3-9, the QPO spends 
less time near the minimum frequency. 
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Fig. 3. — Scatter plots of QPO frequency (in Hz) with power law flux, Flux-PL, (in 
I0~ 8 erg cm _2 s _1 ) for a- and /3-class light-curves. Triangles indicate detections in the first 
100 seconds of entering the dip. Squares are points within the last 100 seconds before exiting 
the dip. The line is a fourth-order best fit to the time evolution of the features. Note that 
for (3-1 through (3-5, a strong linear correlation is apparent. For (3-6, (3-7, and (3-9, the cor- 
relation weakens and we see different degrees of hysteresis. The open circles (non-detections 
of the QPO as defined in Figure 2) do not contribute to the apparent hysteresis pattern. In 
(3-8, the range of frequencies is significantly less, probably due to increased non-detections. 
This is evidence that a continuum of QPO behaviors exists within the /3-class light-curves. 
In a-class light curves, a correlation exists over a smaller range of frequencies as well. 
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Fig. 4. — Scatter plots of QPO frequency (in Hz) with blackbody flux, Flux-BB, (in 
10~ 8 erg cm~ 2 s _1 ) for a- and /3-class light-curves. Symbols are as in Figure 3. A corre- 
lation is observed above 4 Hz, but below 4 Hz, the QPO can change significantly while the 
blackbody flux remains relatively constant. 
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Fig. 5. — Scatter plots of QPO frequency (in Hz) with blackbody temperature, Temp-BB, 
(in keV) for a- and /3-class light-curves. Symbols are as in Figure 3. The a-class light curves 
clearly show no correlation. In most of the /3-class light curves, a correlation is seen but with 
a wider dispersion than that of the power law flux relation. Hysteresis is still apparent in 
(3-7 and (3-9 particularly. 
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Fig. 6. — The time evolution of a Group 1 event. The top panel shows a 1-second resolution 
x-ray light-curve. The next panel shows the fine-binned QPO frequency and the open circles 
are low frequency power peaks where the QPO is not detected. The x-ray power law flux 
(in units of 10~~ 8 erg cm~ 2 s~ 1 ) and blackbody disk temperature are shown at four-second 
resolution. The QPO frequency is most strongly correlated to the power law flux during the 
hard dip, but is also correlated to the blackbody disk temperature. 
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Fig. 7. — The time evolution of a Group 2 event. Panels are as described in Figure 6. 
In this case, the QPO frequency begins rising while the power law flux and the blackbody 
temperature are relatively constant. 
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Fig. 8. — The time evolution of a Group 3 event. Panels are as described in Figure 6. 
The hard dip lasts ~ 1200 seconds — twice the length of the Group 1 and 2 events. The 
QPO frequency behavior is similar in morphology to Group 1, and while strongly correlated 
to the power law features, there is no accompanying correlation to blackbody features. 
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Fig. 9. — The trigger spike of the /3-class light-curves for 1997 August {(3-1 through (3-5) and 
1997 September {(3-6 through /3-9) data at one-second time resolution. The solid lines are 
a double Gaussian plus polynomial fit to the data points. We classify the (3-1 through (3-5 
events as Group 1. These have a strong, narrow, symmetric spike and the underlying flux 
has a positive slope. The (3-6 through (3-9 data we classify as Group 2 events. These have 
weaker, wider, asymmetric spikes and the underlying slope is flat or negative. The (3-8 and 
(3-9 light curves have more disturbed spike morphologies on this time scale. 
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Table 1. XTE Observations 



ID 


XTE DATA ID 


Date Observed 


Start Time (UTC) 


(3-1 


20186-03-03-01 


1997 Aug 14 


04:20:52 


(3-2 


20186-03-03-01 


1997 Aug 14 


05:50:52 


(3-3 


20186-03-03-01 


1997 Aug 14 


07:17:12 


(3-4 


20186-03-03-01 


1997 Aug 14 


09:18:56 


(3-5 


20186-03-03-02 


1997 Aug 15 


07:33:36 


(3-6 


20402-01-45-03 


1997 Sep 09 


06:15:32 


(3-7 


20402-01-45-03 


1997 Sep 09 


08:04:28 


(3-8 


20402-01-45-03 


1997 Sep 09 


09:21:28 


(3-9 


20402-01-45-03 


1997 Sep 09 


09:50:36 


a-10 


50125-01-04-00 


2002 Jul 27 


07:30:50 


a-11 


50125-01-04-00 


2002 Jul 27 


10:24:48 


a-12 


50125-01-05-00 


2002 Jul 28 


06:57:37 



Note. - Observation IDs and dates of the 12 epochs. The 
Start Time is used as a zero-point in the figures unless stated 
otherwise. The IDs are based on the Belloni et al. (2000) classi- 
fications. 
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Table 2. Linear Pearson Correlation Coefficients 



ID 


Group 


TF 


BBN 


BBF 


BBT 


PLN 


PLF 


PLI 


(3-1 


1 


0.96 


-0.62 


0.84 


0.87 


0.87 


0.92 


0.52 


(3-2 


1 


0.93 


-0.63 


0.76 


0.86 


0.87 


0.91 


0.58 


(3-3 


1 


0.96 


-0.62 


0.84 


0.90 


0.81 


0.88 


0.39 


(3-4 


1 


0.96 


-0.66 


0.89 


0.89 


0.88 


0.94 


0.44 


[3-5 


1 


0.96 


-0.57 


0.86 


0.85 


0.90 


0.93 


0.62 


(3-6 


2 


0.90 


-0.53 


0.89 


0.88 


0.79 


0.78 


0.53 


(3-7 


2 


0.76 


-0.42 


0.78 


0.75 


0.63 


0.59 


0.56 


(3-8 


2 


0.84 


-0.40 


0.54 


0.63 


0.68 


0.75 


0.40 


(3-9 


2 


0.87 


-0.36 


0.67 


0.63 


0.76 


0.74 


0.39 


a-10 


3 


0.90 


-0.10 


0.73 


0.48 


0.67 


0.77 


0.47 


a-11 


3 


0.87 


-0.11 


0.57 


0.31 


0.74 


0.83 


0.46 


a-12 


3 


0.87 


-0.11 


0.63 


0.43 


0.64 


0.75 


0.41 



Note. — Correlation of spectral features to the QPO frequency. 
The existence of a correlation is believable for values of \r\ > 0.40 
and highly significant for values of |r| > 0.70. In most cases, the 
correlation to the power law flux is strongest. The abbreviations 
£1X6 clS follows. TF: total flux; BBN: blackbody normalization; 
BBF: blackbody flux; BBT: blackbody temperature; PLN: power 
law normalization; PLF: power law flux; PLI: power law index. 
The Groups are as defined in the text. 
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Table 3. Partial Correlation CoefRcients 



ID 


Group 


PLF — BBF 


PLF BBT 


BBF — PLF 


BBT — PLF 




1 


U.o4 


U.O 1 


U.O 1 


U.O I 


(3-2 


1 


0.84 


0.65 


0.54 


0.37 


(3-3 


1 


0.87 


0.49 


0.82 


0.58 


(3-4 


1 


0.82 


0.74 


0.66 


0.48 


(3-5 


1 


0.84 


0.75 


0.68 


0.39 


(3-6 


2 


0.62 


0.16 


0.81 


0.64 


(3-7 


2 


0.72 


0.08 


0.84 


0.58 


(3-8 


2 


0.80 


0.58 


0.64 


0.31 


(3-9 


2 


0.81 


0.52 


0.77 


0.16 


a-10 


3 


0.80 


0.79 


0.77 


0.54 


a-11 


3 


0.80 


0.83 


0.45 


0.30 


a-12 


3 


0.81 


0.78 


0.72 


0.53 



Note. — Partial correlation of spectral features to the QPO frequency. The 
first column (PLF — BBF) is the partial of the QPO frequency and power law 
flux with the effect of blackbody flux removed. Note that the correlations are 
strong in most of the cases. The second column (PLF — BBT) relates the 
QPO frequency to power law flux, removing blackbody temperature. These 
correlations are weaker, but significant in most Group 1 cases. The third col- 
umn (BBF — PLF) shows the QPO frequency correlation to blackbody flux, 
removing power law flux. Most are believable, suggesting a complex interplay 
between the power law and blackbody features. The fourth column (BBT - 
PLF) shows the QPO frequency correlation to blackbody temperature, remov- 
ing power law flux. In this case, nearly all correlations drop below significance 
showing that power law flux may trace QPO behavior better than blackbody 
temperature. 
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Table 4. Length of Frequency Dip 



ID 


Group 


Min Freq (Hz) 


Time at Min (s) 


Fraction at Min 


(3-1 


1 


2.2 


180 


0.31 


(3-2 


1 


2.2 


240 


0.34 


[3-3 


1 


2.2 


160 


0.24 


(3-4 


1 


2.4 


160(50) 


0.27 


[3-5 


1 


2.2 


180 


0.34 


(3-6 


2 


2.4 


80(20) 


0.15 


(3-7 


2 


2.2 


60 


0.12 


[3-8 


2 


2.2 


100 


0.15 


(3-9 


2 


2.7 


130(0) 


0.21 


a-10 


3 


2.7 


560(0) 


0.50 


a-11 


3 


2.7 


520(0) 


0.44 


a-12 


3 


3.1 


800(0) 


0.63 



Note. - Approximate length of time QPO stays near lowest fre- 
quency. The third column shows the minimum QPO frequency. The 
fourth column shows the length of time the QPO stays within 0.5 Hz 
of the minimum frequency. The number in parentheses is the time 
during which the QPO frequency is below 2.5 Hz if it is not equal to 
the listed time. The fifth column shows the approximate fraction of 
the dip length spent at the minimum frequency. In four cases, the fre- 
quency does not drop below 2.5 Hz. In general, Group 2 events spend 
less time at the minimum frequency than Group 1 events. 
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Table 5. Gaussian Fit Parameters 



ID 


Group 


A 


OA 


Ci 


fint 


(3-1 


1 


0.70 


2.20 


18.24 


82.10 


[3-2 


1 


1.10 


2.39 


9.98 


85.96 


(3-3 


1 


0.90 


2.56 


10.32 


94.75 


(3-4 


1 


1.09 


1.85 


12.14 


89.26 


P-5 


1 


0.86 


2.23 


16.12 


91.85 


(3-6 


2 


0.80 


3.64 


-2.56 


94.49 


(3-7 


2 


0.62 


3.26 


-1.61 


92.97 


[3-8 


2 


0.58 


5.16 


-18.07 


97.10 


(3-9 


2 


0.56 


6.11 


-16.62 


75.23 



Note. - Gaussian fit parameters calcu- 
lated for the "trigger spike" in /5-class light- 
curves. The column labels are as follows: A 
= the normalized amplitude; a a = the width 
(standard deviation) of the gaussian; C\ 
= the slope of the light-curve background; 
fint = the integrated count rate. They are 
grouped by a combination of spike strength, 
spike width, and the sign of the underlying 
slope. 



